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Abstract. - Computer simulations recently revealed that several liquids exhibit strong correlations 
between virial and potential energy equilibrium fluctuations in the NVT ensemble [U. R. Pedersen 
et al, Phys. Rev. Lett. 100, 015701 (2008)]. In order to investigate whether these correlations 
are present also far from equilibrium const ant- volume aging following a temperature down jump 
from equilibrium was simulated for two strongly correlating liquids, an asymmetric dumbbell model 
and Lewis- Wahnstrom OTP, as well as for SPC water that is not strongly correlating. For the two 
strongly correlating liquids virial and potential energy follow each other closely during the aging 
towards equilibrium. For SPC water, on the other hand, virial and potential energy vary with 
little correlation as the system ages towards equilibrium. Further proof that strong pressure-energy 
correlations express a configuration space property comes from monitoring pressure and energy 
during the crystallization (reported here for the first time) of supercooled Lewis-Wahnstrom OTP 
at constant temperature. 



Introduction. — The subject of thermal equilibrium fluctuations is old and well un- 
derstood. For large systems equilibrium fluctuations are approximately Gaussian. For any 
system the "Gaussian components" of the fluctuations - their correlation functions - deter- 
mine the system's response to external fields. This is summarized in the noted fluctuation- 
dissipation theorem known for half a century [1], and little new is expected to be learned 
from studying equilibrium fluctuations 

It was recently shown that several model liquids exhibit strong correlations between 
the thermal equilibrium fluctuations of the configurational parts of pressure and energy at 
constant volume [2-5]. Recall that the pressure p is a sum of the ideal gas term NksT/V 
and a term reflecting the interactions, W/V, where W is the so-called virial [1]: 



pV = Nk B T + W. (1) 

Temperature is defined via the kinetic energy [1] , so the ideal-gas pressure term is a function 
of the particle momenta. The instantaneous virial W is a function of the particle positions, 
W = W(ri, rpf). In the same way, of course, the energy E is a sum of the kinetic energy 
and the potential energy U . If AJ7 is the instantaneous potential energy minus its average 
and AW the same for the virial, at any given state point the WU correlation coefficient R 
is defined by (where sharp brackets denote equilibrium NVT ensemble averages) 
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Strongly correlating liquids by definition have R > 0.9. For simplicity we shall occasionally 
talk about "strong pressure-energy correlations" although it is only the configurational parts 
of pressure and energy - virial and potential energy - that correlate strongly. 

Strongly correlating liquids include [2-5] the standard Lennard- Jones (LJ) liquid (and 
crystal), the Kob- Andersen binary LJ liquid, various binary LJ type mixtures, a dumbbell- 
type liquid of two different L J spheres with fixed bond length [6] , a system with exponential 
repulsion, a seven-site united-atom toluene model, the Lewis- Wahnstrom ortho-terphenyl 
(OTP) model, and an attractive square-well binary model. Liquids that are not strongly 
correlating include water and methanol models [5] . The physical understanding developed in 
Refs. [5,7] is that strong pressure-energy correlations are a property of van der Waals liquids 
and some or all metallic liquids. Liquids with directional bonding like covalcnt and hydrogen- 
bonding liquids do not have strong pressure-energy correlations. Likewise, ionic liquids are 
not expected to be strongly correlating because of the different distance dependence of 
their short-range repulsions and the Coulomb interactions - competing interactions spoil 
the correlations. 

Strongly correlating liquids appear to have simpler physics than liquids in general, an 
observation that has particular significance for the highly viscous phase [8]. It has been 
shown that supercritical (experimental) argon is strongly correlating [2,7], that strongly 
correlating viscous liquids have all eight frequency-dependent thcrmoviscoclastic response 
functions [9, 10] given in terms of just one [3] (are "single-parameter liquids," i.e., have 
dynamic Prigogine-Defay ratio close to unity [7]), that strongly correlating viscous liquids 
obey density scaling, i.e., that their relaxation time r depends on density p and temperature 
as t k F{p 1 /T) [11], and that even complex systems like a biomembranc may exhibit 
significant pressure-energy correlations for their slow degrees of freedom [12]. 

Whenever equilibrium fluctuations of viral are plotted versus those of the potential energy 
for a strongly correlating liquid, an elongated ellipse appears [2,4,5]. The slope 7 of this 
ellipse is given by 



This quantity (that is weakly state-point dependent) is the number entering into the above- 
mentioned density scaling relation [13-15]. This means that for strongly correlating liquids 
knowledge of equilibrium fluctuations at one state point provides a prediction about how 
the relaxation time varies with density and temperature. 

What causes strong WU correlations? A hint comes from the fact that an inverse power- 
law pair potential, v(r) oc r~ n where r is the distance between two particles [16], implies 
perfect WU correlation (R = 1) [2,7]. In this case 7 = n/3. In simulations of the standard 
LJ liquid we found 7 = 6 which corresponds to n = 18. Although this may seem puzzling at 
first sight given the expression defining the LJ potential, Vlj(t) = 4e[(r / 'er)~ 12 — (r/er)~ 6 ], if 
one fits the repulsive part of the LJ potential by an inverse power law, an exponent n = 18 
in fact is required [2,7,17]. This is because the attractive r -6 term makes the repulsion 
steeper than the bare repulsive r~ 12 term would imply. 

Reference [7] gave a thorough discussion of the correlations with a focus on the LJ case, 
including also a treatment of the classical crystal where 0.99 < R < 1 at low temperature 
(an anharmonic effect that survives the T — > limit). According to Ref. [7] the r-dependent 
effective exponent n which controls the correlation is not simply that coming from fitting 
the repulsive part of the potential, but rather n^ 2 \r) = —2 — rv'" (r) / v" (r) . This number is 
approximately 18 around the LJ minimum; in fact the LJ potential may here be fitted very 
well with an "extended" inverse power-law potential [7], vl.t(0 — Ar~ n + B + Cr where 




(3) 
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Fig. 1: Computer simulations of virial and potential energy during the aging of two strongly cor- 
relating liquids following temperature down jumps at constant volume (NVT ensemble) [18]. (a) 
The asymmetric dumbbell model at density p — 1.109 g/ml [6]. The liquid was first equilibrated at 
T=300 K. Here simultaneous values of virial and potential energy are plotted for several times pro- 
ducing the green ellipse, the elongation of which directly shows the strong correlation in equilibrium. 
Temperature was then changed to T—240 K where the red ellipse marks the equilibrium fluctuations 
after equilibration. The aging process itself is given by the blue points. These points follow the 
line defined by the two equilibrium simulations, showing that pressure and energy correlate also 
out of equilibrium, (b) Similar temperature down jump simulation of the Lewis- Wahnstrom OTP 
system [19] with quench data plotted every 0.8 ps. The colors have the same meaning as in (a): 
Green marks the high-temperature equilibrium (T=600 K), red the low-temperature equilibrium 
(T=450 K), and blue the aging towards equilibrium. - Note that in both (a) and (b) the slope of 
the dashed line is not precisely the number 7 of JJq3 (3) ; this is because the liquids are not perfectly 
correlating (the line slope is approximately (AUAW) {(AU) 2 } [5]). (c) Virial and potential energy 
for the asymmetric dumbbell model as functions of time after the temperature jump of (a); in the 
lower subfigure data were averaged over 10 ps. Virial and potential energy clearly correlate closely. 
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Fig. 2: Virial versus potential energy after a temperature down jump at constant volume applied to 
SPC water that is not strongly correlating (colors as in Fig. 1). (a) SPC water at 1 atm equilibrated 
at T=260K, subsequently subjected to a temperature jump down to T=160 K. The data do not 
follow any line, but are pretty much all over the plane, (b) Same simulation at 2GPa where there 
are somewhat stronger WU correlations. The picture is the same, just slightly less confusing. 



n = 18; for this particular potential of course (r) = n. At constant volume the linear 
term contributes little to the viral and potential-energy fluctuations. Thus almost correct 
Boltzmann probability factors are arrived at by using the inverse power-law approximation. 
The equation of state, however, is poorly represented by the inverse power law approximation 
because this potential has no attractive part [7]. 

Temperature down jumps of three model liquids. If the potential at constant 
volume to a good approximation may be replaced by an inverse power-law pair potential, 
strong WU correlations should be present also when the system is not in thermal equilibrium. 
That is, according to the inverse power-law explanation strong WU correlations characterize 
all configurations at a given volume. This can be tested in computer simulations where a 
non-equilibrium situation is easily created by, for instance, abruptly changing temperature 
for a well equilibrated system. 

Figure 1(a) shows the results for a temperature down jump at constant volume, starting 
and ending in equilibrium (NVT ensemble) . The system studied is the asymmetric dumbbell 
liquid consisting of two different sized LJ particles glued together by a bond of fixed length 
[3], a toy model for toluene. The system was first equilibrated at 300 K. The green ellipse 
consists of several simultaneous instantaneous values of U and W in equilibrium at T=300 
K. The strong WU correlation is revealed by the elongation of the ellipse (R = 0.97; 7 = 
6.1). When the liquid is similarly equilibrated at 240 K, the red blob appears. To test for 
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correlation in an out-of-equilibrium situation we changed temperature abruptly from the 
300 K equilibrium to 240 K. The blue points show how virial and potential energy change 
following the temperature down jump. Clearly, strong WU correlations are present also 
during the aging towards equilibrium. Figure 1(b) shows the same phenomenon for the 
Lewis- Wahnstrom (LW) ortho-terphenyl (OTP) model that consists of three Lennard-Jones 
spheres at fixed length and angle with parameters optimized to mimic real OTP [19]. LW 
OTP is also strongly correlating (R = 0.91; 7 = 7.6). The colors are as in Fig. 1(a): Green 
gives the high-temperature equilibrium T = 600 K, red the T = 450 K equilibrium, and the 
blue points show the aging towards equilibrium after changing temperature from 600 K to 
450 K. The picture is the same as in Fig. 1(a): The blue points follow the dashed line. Thus 
virial and potential energy correlate strongly also for far-from-equilibrium states. Figure 
1(c) plots W(t) and U(t) for the Fig. 1(a) data for the asymmetric dumbbell liquid. W(t) 
and U(t) follow each other closely both on picosecond time scales and in their slow, overall 
drift to equilibrium. 

What happens when the same simulation scheme is applied to a liquid that is not strongly 
correlating? An example is SPC water, where the hydrogen bonds are mimicked by Coulomb 
interactions [20] (any model like SPC water with a density maximum cannot have significant 
pressure-energy correlations [2,5]). Figure 2 shows results of simulations of SPC water at 
two different densities, (a) corresponding to low pressure and (b) to very high pressure. In 
the first case virial and potential energy are virtually uncorrelated (R = 0.05 at T=260 
K); in the second case correlations are somewhat stronger (R = 0.34 at T=260 K) though 
still weak. As in Fig. 1 green denotes the initial high-temperature equilibrium, red the 
low-temperature equilibrium, blue the aging towards equilibrium. Clearly, W and U are not 
closely linked to one another. The same is apparent from Fig. 2(b), but here is better WU 
correlation, consistent with the fact that R is larger than in Fig. 2(a). 

Pressure and energy monitored during crystallization of a supercooled liquid: 
The Lewis- Wahnstrom OTP model. — A different far-out-of-equilibrium situation is 
that of crystallization of a supercooled liquid monitored at fixed volume and temperature. 
To the best of our knowledge crystallization of the LW OTP model has not been reported 
before, but Fig. 3 shows that for simulations over microseconds the supercooled liquid crys- 
tallizes at T=375 K and p — 1.135g/cm 3 . The crystal (a distorted fee structure with random 
bond orientations) is shown in Fig. 3(b). Figure 3(a) shows how time-averaged pressure and 
energy develop during crystallization. Contributions to pressure and energy from momenta 
are virtually constant after averaging over 1 ns, so strong WU correlations manifest them- 
selves in strong averaged-pressure /averaged-energy correlations. Clearly averaged pressure 
and averaged energy follow each other closely also during crystallization. This confirms the 
above finding, as well as those of Ref. [5] that strong correlations apply also for the crystalline 
phase of a strongly correlating liquid. Note that the slope 7 is virtually unaffected by the 
crystallization. The persistence of strong pressure-energy correlations during crystallization 
and the insignificant change of the 7's are remarkable, because physical characteristics are 
rarely unaffected by a first-order phase transition. This shows, again, that the property 
of strong pressure-energy correlations pertains to the intermolecular potential, not to the 
particular configurations under study. 

Figure 3(c) shows the radial distribution functions for the liquid and crystalline phases. 

Concluding remarks. — The above out-of-equilibrium simulations show that the 
property of strong WU correlation is not confined to thermal equilibrium. Strongly cor- 
relating liquids have a particularly simple configuration space. These results have signifi- 
cance for any out-of-equilibrium situation. Consider the potential energy landscape picture 
of viscous liquid dynamics [21] according to which each configuration has an underlying 
inherent state defined via a deepest-descent quench, a state that contains most information 
relevant to the slow dynamics (that may be regarded as jumps between different inherent 
states [21]). Figure 4 shows a WU plot of the asymmetric dumbbell model in different 
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Fig. 3: Crystallization of the supercooled Lewis-Wahnstrom ortho-terphenyl (OTP) liquid where 
each molecule consists of three Lennard- Jones spheres with fixed bond lenghts and angles [19]. (a) 
Pressure (right) and energy (left) monitored as functions of time during crystallization at constant 
volume (T= 375 K). Both quantities were averaged over 1 ns; in this way the pressure /energy 
fluctuations directly reflect the virial / potential energy fluctuations. The horizontal dashed lines 
give the liquid (upper line) and crystal (lower line), the averages of which were obtained from 
the simulation by averaging over times 0-2 fis and 5-10 /is, respectively. Both liquid and crystal 
show strong correlations, and the correlations ie^also present during the crystallization, (b) The 
resulting OTP crystal. The crystal is an fcc-like crystal where the LJ atoms form a slightly distorted 
fee lattice with random bond directions, (c) Radial distribution functions of liquid and crystalline 
phases. The two spikes present in both phases come from the fixed bond lenghts. 
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Fig. 4: WU plot for the asymmetric dumbbell model. The upper right corner shows data for 
simultaneous values of virial and potential energy for four equilibrium simulations (T=240-350 K). 
When quenching each of these to zero temperature in order to identify the inherent states, the 
crosses are arrived at. The intermediate points are glasses prepared by different cooling rates, 
i.e., out-of-equilibrium systems generated by cooling in 1 ns from 240 K to the temperature in 
question. This plot shows that strong virial / potential energy correlations are not limited to 
thermal equilibrium situations. 



situations: equilibrium states (upper right) and their corresponding inherent states (lower 
left, one quench per temperature), and glasses at different temperatures in between. A glass 
is an out-of-equilibrium state, and inherent states may be regarded as zero-temperature 
glasses. Altogether this plot shows once again that strong correlations are present also far 
from equilibrium 

If Wis is the inherent state virial, Eq. (1) may be rewritten as follows 

Wis , W-W 1S Nk B T 

p = — + — v— + ^v~- (4) 

The fact that strong pressure-energy correlations apply also out of equilibrium validates 
the argument presented in Ref. [7] relating to the beautiful 2002 paper by Mossa, La Nave, 
Sciortino, and Tartaglia [22], who simulated four different aging scenarious of LW OTP 
(temperature jumps at constant volume, pressure jumps at constant temperature, pressure 
jumps at constant temperature in the glass phase, and isobaric heatings of a glass). Mossa 
et al. showed that all their findings may be rationalized in the following equation of state 
for the averge pressure involving the average inherent energy, eis: 



p(T,V,(e IS )) = Pi S (V,(eis))+p vih (T,V,(e ls )). (5) 

Thus for given temperature and volume knowledge of (eis) in an out-of-equilibrium situation 
is enough to determine the average pressure (we follow the original notation, but eis = 
Uj$). As noted in Ref. [7] Eq. (5) follows from Eq. (4): The pis term derives from the 
Wis/V term of Eq. (4) because (Wis) is a (linear) function of (eis)- For any given state 
W — Wis — l{U — I7is) = l( e — e is) which - assuming thermal equilibrium for the fast 
(vibrational) degrees of freedom - upon averaging implies a contribution to the average 
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pressure that depends on temperature and volume, with only a weak contribution from 
(e IS ). In conjunction with the Nk B T/V term of Eq. (4) this gives the p v ;b term of the 
Mossa et al. out-of-equilibrium equation of state. This derivation, incidentally, explains 
the subsequent finding by Sciortino and coworkers that Eq. (5) does not apply for water 
models [23] - as we have seen water is not strongly correlating. 

In conclusion, strong pressure-energy correlations reside in configuration space. This is 
consistent with the explanation provided in Ref. [7] that the strong correlations reflect the 
fact that an inverse power-law is a good approximation to the potential at fixed-volume con- 
ditions. In particular, the inverse power-law approximation yields almost correct Boltzmann 
probability factors [14]. The inverse power-law approximation is not exact, of course, but it 
reflects a hitherto overlooked approximately scale-invariance of the physics of a large class 
of liquids [14] (van der Waals liquids and, possibly, some or all metallic liquids). Work is 
under way to further investigate the implications of this approximate scale invariance. 

* * * 

The centre for viscous liquid dynamics "Glass and Time" is sponsored by the Danish 
National Research Foundation (DNRF). 

REFERENCES 

[1] L. D. Landau and E. M. Lifshitz, Statistical Physics Part 1 (Pergamon Press, London, 1980); 
J. P. Hansen and I. R McDonald, Theory of Simple Liquids, 2nd ed. (Academic Press, New 
York, 1986); M. P. Allen and D. J. Tildesley, Computer Simulation of Liquids (Oxford Science 
Publications, Oxford, 1993); L. E. Reichl, A Modern Course in Statistical Physics, 2nd ed. 
(Wiley, New York, 1998). 

[2] U. R. Pedersen, N. P. Bailey, T. B. Schr0der, and J. C. Dyre, Phys. Rev. Lett. 100, 015701 
(2008). 

[3] U. R. Pedersen, T. Christensen, T. B. Schr0der, and J. C. Dyre, Phys. Rev. E 77, 011201 
(2008). 

[4] D. Coslovich and C. M. Roland, J. Phys. Chem. B 112, 1329 (2008). 

[5] N. P. Bailey, U. R. Pedersen, N. Gnan, T. B. Schr0der, and J. C. Dyre, J. Chem. Phys. 129, 

184507 (2008). 

[6] A system consisting of 512 asymmetric dumbbell molecules modelled as two Lennard- Jones 
spheres connected by a rigid bond was simulated. The dumbbells were parameterized to mimic 
toluene. A large sphere (mimicking a phenyl group) was taken from the Lewis- Wahnstrom OTP 
model [19] with the parameters m p — 77.106 u, a p = 0.4963 nm and e p — 5.726 kj/mol. A 
small sphere (mimicking a methyl group) was taken from UA-OPLS having m m = 15.035 u, 
o~ m = 0.3910 nm and e m = 0.66944 kj/mol. The bonds were kept rigid with a bond length of 
d = 0.29 nm. The volume was V = 77.27 nm giving an average pressure of approximately 1 
atm. The temperature was held constant at T = 130 K using the Nose-Hoover thermostat. 

[7] N. P. Bailey, U. R. Pedersen, N. Gnan, T. B. Schr0der, and J. C. Dyre, J. Chem. Phys. 129, 

184508 (2008). 

[8] W. Kauzmann, Chem. Rev. 43, 219 (1948); G. Harrison, The Dynamic Properties of Su- 
percooled Liquids (Academic, New York, 1976); S. Brawer, Relaxation in viscous liquids and 
glasses (American Ceramic Society, Columbus, OH, 1985); I. Gutzow and J. Schmelzer, The 
Vitreous State: Thermodynamics, Structure, Rheology, and Crystallization (Springer, Berlin, 
1995); M. D. Ediger, C. A. Angell, and S. R. Nagel, J. Phys. Chem. 100, 13200 (1996); R. V. 
Chamberlin, Phase Transitions 65, 169 (1999); C. A. Angell, K. L. Ngai, G. B. McKenna, P. 
F. McMillan, and S. W. Martin, J. Appl. Phys. 88, 3113 (2000); C. Alba-Simionesco, C. R. 
Acad. Sci. Paris (Ser. IV) 2, 203 (2001); P G. Debenedetti and F. H. Stillinger, Nature 410, 
259 (2001); K. Binder and W. Kob, Glassy Materials and Disordered Solids: An Introduction 
to their Statistical Mechanics (World Scientific, Singapore, 2005); F. Sciortino, J. Stat. Mech., 
P05015 (2005); J. C. Dyre, Rev. Mod. Phys. 78, 953 (2006). 

[9] N. L. Ellegaard, T. Christensen, P. V. Christiansen, N. B. Olsen, U. R. Pedersen, T. B. 
Schr0der, and J. C. Dyre, J. Chem. Phys. 126, 074502 (2007); N. P. Bailey, T. Christensen, 



p-8 



Strong pressure-energy correlations as a configuration space property 



B. Jakobsen, K. Niss, N. B. Olsen, U. R. Pedersen, T. B. Schr0der, and J. C. Dyre, J. Phys.: 
Condens. Matter 20, 244113 (2008). 
[10] T. Christensen and J. C. Dyre, Phys. Rev. E 78, 021501 (2008). 

[11] A. Tolle, Rep. Prog. Phys. 64, 1473 (2001); C. Dreyfus, A. Aouadi, J. Gapinski, M. Matos- 
Lopes, W. Steffen, A. Patkowski, R. M. Pick, Phys. Rev. E 68, 011204 (2003); C. Alba- 
Simionesco, A. Cailliaux, A. Alegria, and G. Tarjus, Europhys. Lett. 68, 58 (2004); R. Casalini 
and C. M. Roland, Phys. Rev. E 69, 062501 (2004); C. M. Roland, S. Hensel-Bielowka, M. 
Paluch, and R. Casalini, Rep. Prog. Phys. 68, 1405 (2005). 

[12] U. R. Pedersen, G. H. Peters, T. B. Schr0der, and J. C. Dyre, arXiv:081 1.3317 (2008). 

[13] T. B. Schr0der, U. R. Pedersen, and J. C. Dyre, arXiv:0803.2199 (2008). 

[14] T. B. Schr0der, U. R. Pedersen, N. P. Bailey, S. Toxvaerd, and J. C. Dyre, arXiv:0812.4960 
(2008). 

[15] D. Coslovich and C. M. Roland, J. Chem. Phys. 130, 014508 (2009). 

[16] W. G. Hoover and M. Ross, Contemp. Phys. 12, 339 (1971); W. G. Hoover, D. A. Young, and 
E Grover, J. Chem. Phys. 56, 2207 (1972); Y. Hiwatari, H. Matsuda, T. Ogawa, N. Ogita, 
and A. Ueda, Prog. Theor. Phys. 52, 1105 (1974); L. V. Woodcock, Phys. Rev. Lett. 54, 1513 
(1985); P. G. Debenedetti, F. H. Stillinger, T. M. Truskett, and C. J. Roberts, J. Phys. Chem. 
103, 7390 (1999); E. La Nave, F. Sciortino, P. Tartaglia, M. S. Shell, and P. G. Debenedetti, 
Phys. Rev. E 68, 032103 (2003); G. Rickayzen and D. M. Heyes, Phys. Rev. E 71, 061204 
(2005); R. Casalini, U. Mohanty, and C. M. Roland, J. Chem. Phys. 125, 014505 (2006). 

[17] D. Ben-Amotz and G. J. Stell, J. Chem. Phys. 119, 10777 (2003). 

[18] NVT simulations were carried out using Gromacs software [H. J. C. Berendsen, D, van der 
Spoel, and R. van Drunen, Comp. Phys. Comm. 91, 43 (1995); E. Lindahl, B. Hess, and D. 
van der Spoel, J. Mol. Mod. 7, 306 (2001)] using the Nose-Hoover thermostat [S. A. Nose, 
Mol. Phys. 52, 255 (1984); W. G. Hoover, Phys. Rev. A 31, 1695 (1985)]. Molecules were kept 
rigid using the LINCS algorithm [B. Hess, H. Bekker, H. J. C. Berendsen and J. G. E. M. 
Fraaije, J. Comp. Chem. 18, 1463 (1997)]. 

[19] L. J. Lewis and G. Wahnstrom, Phys. Rev. E 50, 3865 (1994). 

[20] H. J. C. Berendsen, J. R. Grigera, and T. P. Straatsma, J. Phys. Chem. 91, 6269 (1987). 
[21] M. Goldstein, J. Chem. Phys. 51, 3728 (1969); F. H. Stillinger and T. A. Weber, Phys. Rev. 

A 28, 2408 (1983); F. H. Stillinger, Science 267, 1935 (1995); T. B. Schr0dcr, S. Sastry, J. C. 

Dyre, and S. C. Glotzer, J. Chem. Phys. 112, 9834 (2000). 
[22] S. Mossa, E. La Nave, F. Sciortino, and P. Tartaglia, Eur. Phys. J. B 30, 351 (2002). 
[23] N. Giovambattista, H. E. Stanley, and F. Sciortino, Phys. Rev. Lett. 91, 115504 (2003). 



p-9 



